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A central issue in the ongoing studies of graphene is how the morphology of the sample 
modifies its electronic properties. For curved shapes its well-known special relativistic-like 
properties naturally point to quantum field theory in curved spaces that is not a general 
relativistic-like description as yet because the time components of the metric are left over. 
Here we find that, if graphene is shaped as Beltrami pseudosphere the problem is solved. As a 
result, through the appropriate identifications of the graphene counterpart of the customary 
relativistic procedures, we are able to predict that the curved graphene sample has a finite 
temperature electronic local density of states that is a realization of the Hawking-Unruh 
effect. We then propose an experiment with a Scanning Tunneling Microscope that should 
easily detect such effect and measure the Hawking temperature. 

1 Graphene and quantum field theory in curved spacetimes 

Graphene is an allotrope of carbon that was first theoretically positecP^ and then found to have 
an abundance of "unorthodox" propertieP'^, the understanding of which is appealing to theorists 



and experimentalists, see, e.g.,^. Following the spirit of the condensed matter analogues of grav- 
itational effectPI^ we shall use graphene as a physical realization of quantum field theory (QFT) 
in curved spacetimeP^H] Qf these unorthodox properties one that is crucial is the dimensionality 
(we shall evoke the exotic scenarios of gravity in 2+1 dimension^^^Hl) and another one is the well 
known massless Dirac field description of graphene's conductivity electrons^"^ (we shall use the 
fact that the theory enjoys a scale symmetry proposed by Weyl in an early attempt to introduce the 
gauge principlePE^. 

Li this work we assume that the conductivity electrons of graphene experience a general 
relativistic-like spacetime and we identify a specific shape, the Beltrami pseudosphere, for which 
it is easier to probe this assumption. In turn, the predictions of QFT in curved spacetimes point to 
a thermal spectrum in the form of a finite temperature electronic local density of states (LDOS). 
The temperature is of the Hawking-Unruh type ElES of geometric origin and depends upon 
the curvature and upon the the meridian coordinate on the surface. The final outcome of this 
study is an exact formula for the LDOS that we show to be possible to measure. The whole 
construction, including the assumption on graphene's electrons experiencing a curved spacetime, 
can be submitted to the judgement of experiments. Among the unusual features of the formula is a 
Bose-Einstein thermal spectrum for this Dirac system that is due to the odd dimension^^SIIMQ] 
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2 Shaping graphene as a Beltrami pseudosphere and the issue of time 



For graphene, to include time in a relativistic fashion we need to take two steps. First, we have 
to assume that for a curved graphene sample the conductivity electrons experience a spacetime 



metric, the choice dictated by everyday practice being 



graphene 



1 




dap 



V 







(1) 



i.e., the time-time component is just the flat one. Here z/ = 0, 1, 2 and a, (3 = 1,2, = (t, u, v) 
where t coincides with the laboratory time and u, v are the coordinates on the surface. Second, we 
have to use a Lagrangian rather than a Hamiltonian description. Consequently, the dynamics of 
graphene's conductivity electrons near a Dirac point is given by the customary generalization to a 
curved spacetimeP^of the action for massless Dirac spinors in 2+1 dimensiond^ (Supplementary 



Information I) 



(2) 



where h = vp = = I, with vp the Fermi velocity. 



We are building upon the continuum description of a Dirac quantum field, and upon modeling 
the effects of curvature through the coupling of the Dirac field to a curved spatial metric. The 
effectiveness of both assumptions was proved elsewhereP^. Here we assert that the electrons 
on graphene might directly experience a curved spacerime even though the curvature is all in the 
spatial part. That a spacetime can be curved acting solely on the spatial part is well knworP', the 



breakthrough is to prove that this makes sense for graphene. In what follows we shall extract a 
measurable prediction based on such assumption and we shall show that the physical set-up to do 
that boils down to shape the graphene sample as a Beltrami pseudosphereP^. 

When the metric ([T]) is conformally flat (Supplementary Information I), we can make use of 
the Weyl symmetr}^^ enjoyed by the action Q to obtain exact results that otherwise are difficult 
or impossible to obtain. As proved in earlier work'^ (see also Supplementary Information I), 
the metric ([T]) is conformally flat for all surfaces of constant Gaussian curvature /C, but we need 
(HI) to explicitly take the form 5'^IfP'^'''''=(Q) = ^'^{Q)glf{Q), for certain coordinates Q^. Only 
then we can use all the power of Weyl symmetry. We know from general arguments that those 
coordinates must exist, but the key issue for graphene is whether such coordinates Q'^ can be 
practically realized in the laboratory. 

One of the main results of this work is that the Beltrami pseudosphereP^, M'^ = gai3(f(l^ = 
du^ + r^e^^/^ciw^ with v E [0, 27r], u G [—00, 0] (see Fig. 1), solves the problem. For this surface 
the coordinates are the coordinates q^^, and the time coordinate coincides with the laboratory 
time. The isothermal coordinates are x = v/r and y = e~"/'^ /r (Supplementary Information II), 
hence 

^^raphenc = dsfs) = e'^^^ [e-'^^^dt^ - du^) - T^dv^] . (3) 

As the line element in square brackets is flat (and Rindlei'^'^, an exciting instance that we shall 
soon exploit), for a Beltrami pseudosphere we can fulfill the condition of a physically doable 
gf^^^^'^^{q) = ^'^{(i)g^{(l) already in the frame q^. From now on, the spacetime we shall suppose 
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to be experienced by graphene's quasi-particles, is 



1 







9, 




-1 







(4) 
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where, t E [—oo, +oo], u E [—oo, 0], v E [0, 27r]. We call this a "Beltrami spacetime". 

Once we take this view, any measurement is an operation relating different spacetimes: the 
inner (2+l)-dimensional curved spacetime and the outer (3+l)-dimensional flat spacetime. See 
Fig. 2. The electrons, far from the graphene surface, loose the lattice-induced properties of quasi- 
particles, i.e. electric neutrality, masslessness, (2-1-1 )-dimensionality and, the pseudo-relativistic 
Dirac nature. This latter fact, in particular, requires the final ground state reached after the mea- 
surement be of a non-relativistic nature. To accommodate this hybrid situation into the fully rela- 
tivistic scenarios of QFT in curved spacetimes, when we take into account the different quantum 
vacua for the different observers involvecP'^HnHHESl^ we approximate the ground state associated 
with the measurements with the Minkowskian \0m)- This way we have both, the mathematical 
description of the experimental evidences of the lattice-induced features (we keep a relativistic- 
like, (2-i-l)-dimensional structure), and the information on the true ground state (we use a flat 
vacuum defined everywhere). This is the best approximation: the Dirac field is living in a (2-1-1)- 
dimensional curved spacetime with coordinates = {t, u, v), the measuring device has the same 
coordinates (i.e. it follows the profile of the surface in a specific manner described later), the quan- 
tum vacuum of reference entails information on the flatness of the ambient spacetime but retains 
the basic lattice-induced features. 
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3 Weyl symmetry and the fictitious Rindler spacetime 



We shall focus on the one particle Green's function that contains all the information on the single 



particle properties of the system such as the LDOS, life time of the quasi-particles and thermody- 



namic properties (specific heat). For the reasons illustrated above, this is defined as 



(5) 



i.e. the positive frequency Wightman function, in the language of QFT on curved spacetimes^'^^ 
(see also^^. To obtain an exact result for this function, we shall use local Weyl symmetrj^^G^ as 
this case is a perfect match for its implementation: (7^^^ = ip'^{u)g^^\ ip^^^ = ip~^{u)4j^^\ with 
ip{u) = e"/^ and the metric 







describing a flat geometry. Local Weyl symmetry gives 
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Notice that Ab refers to a curved spacetime, whereas Ar to flat curvilinear coordinates. The 



Green's function of interest can then be written as 



(8) 



where 



(9) 



The superscript R here stands for "Rindler". We shall now study the spacetime gliu to compute 
5**^^^ but one must bear in mind that this is a fictitious flat spacetime, the only physical spacetime 
being the (Weyl-equivalent) curved Beltrami spacetime gl^\ The Weyl transformation connects 
the Beltrami spacetime to a flat spacetime that, as well known, is the near-horizon spacetime for a 
black hole--. We claim that this setting is the closest possible to a black hole on graphene that can 
be obtained by solely acting on the spatial part of the metric. The physical Green's function S^^^ 
is recovered from S^^^ through ([8]). 

Let us introduce Minkowski coordinates = (T, X, Y) for which (is^^^ = g\^dq^dq'^ = 
fl^vdQ^dQ": T = re~"/^sinh^, X = rv, Y = re""/'' cosh ^. Although no physical quantity 
explicitly depends from these coordinates, we have assumed that the quantum vacuum of reference 
is Minkowskian. We have that — = r^e~^"/^' = a~'^{u), which, for constant u, corresponds 
to worldlines of observers moving at constant "proper acceleration' SUOHI] 

a{u) = ^ E[0,r-^]. (10) 

In a standard Rindler spacetime, the acceleration a can become infinite. In our case there is a 
maximal "acceleration", amax = a{u = 0) = = \J —K, of complete geometric origin. Fur- 
thermore, we are forever confined to one Rindler wedge (a > 0) and on one side (n < 0) of the 
singular maximal circle ("Hilbert horizon"), the other side is unaccessible as "the world ends" at 
R = r which corresponds to n = (see z{u) in the caption of Fig. 1). For a radius of curvature 
r ~ 1mm, the Rindler horizon F = T is effectively reached after a time thor ~ r/vp of few 
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nanoseconds (vp ~ 10 m/s). Let us also introduce a perhaps more familiar notation 

7] = - = a^^J, ^ = re""/^ = , (11) 

r 

with which c/s^^^ = ^'^drj'^ — d^'^ — r'^dv'^. 

The Rindler spacetime has been studied in great depth, many crucial results for massless 
fields are exacP^, and is well known to give rise to a thermal spectrum for quantum field^. Our 
task is to see whether those results can apply to graphene. The real coordinates g'^ are, through (fTTl) . 
those of an "accelerated" observer in the fictitious spacetime and constant "acceleration" means 
u = constant = u (we also take v = constant = v). Thus, measuring at a fixed point (n, v) 
of the real Beltrami surface and letting time pass corresponds, in the equivalent fictitious Rindler 
spacetime, to measurements made by an observer on a worldline of constant proper "acceleration" 
a: t/r = a{u)T, u = u,v = v. 

To fit our set-up within the proper relativistic requirements we need then to refer to the 
Green's function ([9]) at the same point in space and at two dijferent times 

S^^'^t, - t2, q, q) = (OmI^^^H^i; ^)^^''\t2; q)|OM) , (12) 

as this is what would be seen by an observer of the above described worldline. The dependence on 
ti — t2 is a result of the stationarity of the worldline in point. Furthermore, we need to compute 
Si^) by setting t ^ t + is (where t\ — t2 = t) di prescription that takes into account the nonzero 
size of the measuring deviceP^ni. For the Scanning Tunneling Microscopy (STM) experiment we 
have in mind, e is the size, in "natural units", of its needle or tip. For a tungsten needle e ~ 
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0.25mm x Vp^ ~ 10~^°s, while for a typical tip e ~ lOA x Vp^ ~ 10~^^s (see, e.g.j^. Hence, the 
Unruh requests are satisfied by considering the Green's function S^^^r, q, q), where r = t/e"/'', 
and measuring at each point for a time given by the largest among e and the thor (that is ~ lO^^s 
for r ~ 10mm and ~ 10~^^s for r ~ Ijim). The Unruh thermal features are then readily seen by 
considering, as customary, the power spectrurrpESl 

r+oo 

(13) 



1 r I'+oo 
F(^) (u, q) ^ -Tr 7° J dre-'^^S^^^ (r, q, q) 



that, for graphene, besides inessential constants, coincides with the definition of the electronic 
p'^^\(jj, q) = -F^^\u, q). As we are in a massless model, the quantity in (fT3l) can 
be computed exactly to giveP'F(^)(a;, q) = |c<;/(e'^/^ — 1), where T is a Unruh temperature 

r.^ = ^^.r„e«/', (14) 

that includes the constant % = 1/ (27rr) and the Tolman factoi'^ e^^^, required by the local mea- 
surements. We see here the appearance of a thermal distribution, the typical effect predicted by 
QFT in this context. That is a Bose-Einstein distribution for a Dirac field, an instance theoretically 
predicted to be a general feature ("swapping of statistics") of both the Hawking and Unruh phe- 
nomena in odd dimensional spacetime^^^Clllini Nonetheless, as said in various occasions, the only 
physical set-up we have is the Beltrami spacetime. To discuss the physics we have now to move to 
that case. 
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4 Experimental prediction of the physical LDOS 



The task is easy due to Weyl symmetry. In fact, we have tuned from the very beginning the 
whole procedure towards the optimal use of this symmetry. The Weyl factor in ([8]) being time- 
independent, goes through the Fourier transform in (fT3]) . hence the LDOS we predict for a graphene 
sheet shaped as a Beltrami pseudosphere is p^^^co, q) = (p~^(q)p*^^^(a;, q), or, more explicitly 

4 1 E g"^^/'' 

pW(E, u, r) = -— — , (15) 

71 [nvF) exp [£//(fcB7oe"/'^)J - 1 

where we included the g = A degeneracy, and the proper dimensional units are re-introduced, for 

instance uj = co/vp, E = Jtlo and % = hvp / {kBl-nr). This equation is our prediction of the 

exact functional dependence of the LDOS from the energy E, the meridian coordinate u and the 

curvature r (even if we explicitly wrote-in To this is essentially ~ r~^). With the notation u we 

want to stress that: i) the measuring device has to closely follow the profile of the surface, so that 

its coordinates can be taken to be q^; and ii) the device has to stop at each given point on the surface 

(m, v) for a time much bigger then max(£:, thor). not a stringent condition, as we saw earlier. 

We shall soon describe how to detect this Hawking-Unruh effect through an STM experi- 
ment, but let us first make some comments. We do not have a black- hole here, but we do not need 
one as our Hawking effect is inferred from an equivalent (through Weyl symmetry) Unruh phe- 
nomenon. Similar in spirit are the well known results of where the Hawking effect is obtained 
from a higher dimensional Unruh effect. Our results differ from those because we ruled-out the 
embedding into un-physical higher dimensional spacetimes and because the Beltrami spacetime is 
not a proper anti de Sitter (AdS) (as a consequence, our results are closer to the dS case of rather 
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thanks AdS). 



Let us focus now on the experimental detection. To realize with an STM device in the lab 
the situation depicted in Fig. 2 we need to follow the prescriptions i) and ii) above and we have to 
make sure that the tunneling current measured is the one relative to the sample electrons tunneling 
to the tip, and not viceversa. This is easily achieved by fixing the right polarity of the bias voltage. 
The plot we need is a current-to-voltage (dl/dV) spectrum. This is obtained by first constructing a 
one-layer graphene sample shaped as a given Beltrami pseudosphere with fixed r and by measuring 
at a particular fixed value of u. We then vary the bias voltage to obtain the p^^^ vs E behavior. As 
the formula (fT5l) does not depend on the angle v, we do not need to stay at a fixed point (m, v) on 
the surface but we can loop around the circle at fixed u (see Fig. 1) while varying the bias voltage 
V. We can then repeat the same measurements at different values of u to obtain a series of dl/dV 
spectra (a "line-cut"). The expected line-cut is illustrated in Fig. 3 for r = 10/im. 

To avoid defects proliferation we need to take large curvature on the natural scale furnished 
by the lattice spacing i ~A : r >> £. We take r ~ 10~^m as the highest curvature. This 
corresponds to To ~ 13K as the highest detectable temperature (this value can only be made 
smaller by including the effects of the local measurement, e"/'', as n < 0). The zero-curvature 
limit (large r) of p^^^^ does not match the flat LDOS, = This is as it must be 

since (fTS)) is not the result of a perturbative computation with p*^^^'^ (E) as the leading term, thus 
even a very small curvature fully turns on the effect. In Fig. 4 we show the behavior of p^^^ for a 
wide range of Beltrami pseudospheres with curvatures spanning a range of 6 orders of magnitude. 
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As can be seen there, fixing the value of u in units of r, the effect is the same for all curvatures. 

The limit for zero energy gives a nonzero result, p^^^ (0, u, r) = as expected for a 

finite temperature LDOS. A direct measurement of this value of the LDOS would give the constant 
Hawking temperature % = f ^e"/^p'^^^(0, u, r), nonetheless, this measurement seems difficult to 
perform due to the extremely high sensitivity necessary for the appreciation of this effect, see inset 
of Fig. 4. It would be difficult to realize whether it is the effect of a finite temperature or, e.g., of 
impurities. 

The best way to measure the Hawking temperature is the one based on the line-cut in Fig. 3 
and described above. For each power spectrum (fixed u) we shall obtain the corresponding u- 
dependent temperature, T = -E/ [fcg In ( ^^^^^ ^^pW + 1)]. Different Ts for different us will then 
give the same constant Hawking temperature through % = e~"/^ T. 
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SUPPLEMENTARY INFORMATION 



The Hawking-Unruh phenomenon on graphene 

A. lorio, G. Lambiase 

I. Graphene Dirac spinors and geometric conventions. Due to the special geometry of 
graphene's Honeycomb lattice (two interpenetrating triangular lattices), the Bloch wave function 
has two components representing the probability amplitude at each sublattice [SI]. The spinorial 
nature of the excitations is linked to this sublattice degree of freedom named pseudospin in the 
graphene literature. The other important ingredient of the Dirac model, also a consequence of the 
special topology of the lattice, is the fact that the tight binding dispersion relation at half filling has 
six Fermi points instead of a Fermi line [SI]. As only two of them are inequivalent this introduces 
an extra degeneracy of the vacuum. This new degree of freedom is called valley in semiconductor 
language and plays an important role in the physics of graphene [SI]. Linearizing around the two 
inequivalent Fermi points (Dirac points), k^, k± ^ + p, the tight-binding Hamiltonian can be 
written as (/i = 1) [SI] 



with a = ((7i,cr2), a* = (— 0-1,0-2), o-j being the Pauli matrices, the Fermi velocity vp = 3ti/2 
(that will be set to 1) t ^ 2.7 eV is the hopping parameter, £ ~ 2.5A is the lattice spacing. 
Where ip"^ = (a+ 6+), ip"^ = (a_ 6_) are two-component Dirac spinors, as appropriate for this 



H 




p 




(SI) 
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2+1 -dimensional system, and a and b are anti-commuting annihilation operators for an electron 
in the two sub-lattices. In here, we do not consider short range scattering centers or any other 
effect mixing the two Fermi points, thus we discuss the physics around a single Fermi point, e.g. 

tp = tjj^. The corresponding action is A = i J (Px^'y°'da i), where 7° = o"3, 7^ = ia2, 7^ = —iai 
which obey [7*^, 7*]+ = 2ri'^^, with a, 6 = 0, 1, 2 the Lorentz/flat indices (see below). The curved 
action is 



e^-E^ = 5^, the indices = 0, 1, 2 respond to diffeomorphisms (Einstein indices), while a,b = 
0, 1, 2 respond to flat space transformations (Lorentz indices), rjab = diag(+l, —1,-1), = 
det g^^y and the diffeomorphic covariant derivative is = (S^ + ^w^^'^Jtc), with J"'' = ;|[7",7^], 
and ilJ^^\ = e'{{5^d^ + r^^)i?^ is the spin connection obtained by requiring the full covariant 
derivative of the Vielbein to be zero (metricity condition) V^e" = d^el — F^^e" + ^e]^ = 
0, where F^^^ is the Christoffel connection. Our convention for the Riemann curvature tensor is 
R^xixv = d[yTPf,]x + F^^^^F'^^iA, that, in terms of u^l is {R^uTh = {d[u to,,] + uj[u (^1,])"^ Torsion 
T^^j^ = T'^j^j^j is taken to be zero. 

For the metric (H]) the Ricci tensor is = diag(0, /C, /C) which gives as the only nonzero 
components of the Cotton tensor, C"" = e'^'PVaRp'' + H ^ i^, the result C^^ = -dyJC = and 
C^y = dxJC = C^°. Since conformal flatness in (2+l)-dimensions amounts to C^'^ = 0, this shows 
that all surfaces of constant JC give rise in ([T]) to conformally flat 2-1-1 dimensional spacetimes. 
Note that the result holds for (+, — , — ) and for (+, +, +). 




(S2) 



where we implicitly introduced the Vielbein and its inverse E^, riahe^^e\, = g,,^, e' 
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II. Beltrami spacetime and physical coordinates Q^. We found that the coordinates 
for the sphere, the first example one would consider, are not easy to envisage, while surfaces 
with constant negative curvature, /C < 0, are a better candidate. For the latter case the metric of 
graphene can be written in the spatial isothermal coordinates [S2], i.e., in terms of the line element, 

as (iSgj,3^pjjgjjg = dt^ — ^{dx^ + dy"^) where x, y are the abstract coordinates of the Lobachevsky 



geometry in the upper half -plane (y > 0) model and r = a/— [S3]. The actual type of surface 
is obtained by specifying x and y in terms of coordinates measurable using the Euclidean distance 
(embedding). This is crucial for the metric to describe the spacetime experienced by the electrons 
on real graphene (an instance not considered in [S4,S5]). 

The key fact common to all these surfaces (whether or not they are surfaces of revolution) is 
a theorem by Hilbert that proves the impossibility to have a complete surface of constant negative 
Gaussian curvature embedded in [S6]. Complete means, essentially, that there are no singu- 
larities (if one changes the dimensions, e.g., — M.*, or the signature, (+ + +) — ( — h +), 
the task becomes indeed possible), hence we must expect singularities of some sort in our sur- 



Kdf - dx^ - dy' 



it is easy to check 



faces. By writing the line element as ds^^^^^^^^^ = ^ 
that the line element in square brackets is flat. This apparently solves our problem as the 
that should make sense appear to be {t, x, y). But it is not so until we identify the specific sur- 
face and reduce to the measurable coordinates. For instance, we found that the elliptic pseu- 
dosphere [S6], df^ = du^ + r"^ siwh^ u/rdv"^ (here u and v are the meridian and parallel co- 
ordinate on the surface, respectively) has as natural coordinate system a different one, namely 
(T = e*/'' coshw/r, X = e*/^ sinhw/r cosw, Y = e*/^ sinhw/r sinw ). Clearly, it is very difficult to 
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imagine how to practically realize a physical situation in the laboratory so that those coordinates 
become meaningful. It should be now evident that the identification of the right surface is quite an 
art. 
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Figure 1 Beltrami pseudosphere. The coordinates of the Beltrami pseudosphere, 

in the canonical form, {R{u) cosv, R{u) sin v,z{u)), are x{u,v) = ce"/''cosf, y{u,v) = 
ce"/''sinw, z{u) = r{^l- fle^^/'^-arctanhy^l - ^e^^A), c> and r = ^/-}C~^ > with /C 
the Gaussian curvature. We choose c = r, thus R{u) = re"/*" g [0,r] as n g [-cx),0]. The 
curves at constant angle v are the meridian curves and coincide with the null geodesies 
of the spacetime {d£^ = dt^). 

Figure 2 Choice of vacuum. Graphene quasi-particles enter the measuring device as 
a (2+1)-dimensional massless Dirac "curved" field ^^^^ (a). Eventually, they loose the 
lattice-induced features, and behave as charged, massive, non-relativistic electrons in a 
(3+1)-dimensional spacetime (c). To include both, the routinely detected lattice-induced 
features, and the flat status of the ground state eventually reached, the Green's functions 
need to be expectations of products of ^^^'^ fields on a flat (2-i-1)-dimensional vacuum 
|0m) (b). The IVlinkowski vacuum, being defined everywhere, mimics the true vacuum \Qe) 
better than the Rindler vacuum |0/j) that only lives in one sector of the spacetime. 

Figure 3 Line cut. Series of dl/dV spectra ("line-cut") obtainable with an STIVI. In this 
example the curvature is r = lO/xm. The measurements need to be performed at a fixed 
value of u, but since p^^^ does not depend on the angle v, we can loop around the circle 
at fixed u while varying the bias voltage V. We can then repeat the same measurements 
at different values of u to obtain the shown line-cut. The polarity of the bias V is such 
that the sample electrons tunnel to the STIVI tip, and not viceversa (see Fig. 2). The 
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energy range \s E e [-1, l]eV, where the linear behavior of the flat LDOS (shown in blue) 
is more pronounced. The range of the meridian coordinate is u e]0,-3r]. The u = 
value (red dashed curve) corresponds to the singular boundary of the pseudosphere, 
hence we shall not trust our model there. The u = -3r value (green curve) is chosen 
because beyond that value the pseudosphere is too sharp (line shaped, see Fig. 1). To 
each spectrum corresponds a w-dependent temperature T that ranges from ~ 0.05mK till 
~ O.IK. The constant temperature (same for the whole line-cut) is then obtained through 
% = e-"/'^r ~ 0.13K. 



Figure 4 Nonperturbative effect. dl/dV spectra obtainable with an STIVI for different 
Beltrami pseudospheres (different rs). The main plot shows that the effect is present for a 
wide range of curvatures, from very large to very small r e [10"'^, 10~^]m. Here the value 
of the meridian coordinate is chosen to be u = -r for all curvatures, giving rise to the 
same curve for them all. The inset shows the near zero behavior of p^^^ for the highest 
curvature considered here, r = lO^'^m, that is the most favorable for the experimental 
detection of a finite temperature at E = 0. The intercept of dashed lines with the y-ax\s 
gives the p^^'^{E = 0) values. Despite of that the required instrument sensitivity for such 
detection is too high (the axes of the inset have the same units as for the main plot), and 
it would not be clear whether it is the effect of a finite temperature or, e.g., of impurities. 
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